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1. Introduction 


The computer code BDSCAN, written in FORTRAN by Bruce H. Faaland to 
implement the bound-and-scan algorithm developed by Frederick S. Hillier [1], 
is documented and listed in this report. This program is designed to solve the 


pure integer linear programming problem, 


n 
maximize Xo = C.x, 
j=l | 
subject to 
n 
(1) a ay 5%; <4 (4,52... fu); 
j=1 
(2) Bs a 0 (j = WAP cee) a 
(3) *, is an integer Gira elec Gn) 


This problem may also be written as maximize cx subject to Ax < b, x > 0, 
x integer where the constraint matrix A is mxn, the cost coefficients c 
form a 1xn_ row vector, the right hand side b is an mxJ1 column vector 
and x is an n x1 _ column vector. 

It is required that a suboptimal feasible (integer) solution has been 
obtained by a suboptimal algorithm such as that developed by Hillier [3] and 
documented in a companion report. See [2] for a detailed discussion of the 


bound-and-scan algorithm. 
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2. Main Program and Input Requirements 

The main program for the bound~-and-scan code BDSCAN serves to coordinate 
the 16 parts described in detail in [2] and summarized in [1]. The first card 
of input must contain NAME, any alphanumeric characters to identify the problem. 


The second card contains M, N and KL in 315 format. These variables have the 


following significance: 


M = the number of rows in A. 


N = the number of columns in A. 
1, if the optimal basis inverse is in normalized form 


0, otherwise, 


This basis inverse corresponds to an optimal solution to the stated problem 
ignoring the integer requirements (the "LP problem"). It is in normalized 
form if the constraints Ax < b were normalized as defined in [3, p. 605] 
before the LP problem was solved. 

The third group of cards contains the arrays A, b and c in 15F5.0 
format. The A matrix is read in one row at a time, then b is read in 
Starting on a new card, followed by c also starting on a new card. 

The fourth group of cards contains the basis inverse for an optimal 
solution to the LP problem; the format is 6F13.5. This matrix is read in row 
by row. 

The fifth group of cards contains JPM(I), I = 1,..,M in 1514 format where 
JPM(1) is the index of the ee basic variable (including slacks) in the 
solution to the LP problem. 

The sixth group of cards contains an optimal solution to the LP problem in 
6F13.5 format. The final group of cards contains a suboptimal feasible solution 
to the ILP problem (the LP problem plus the integer requirements) in 6F13.5 
format. 

Another problem may then be input using the card sequence described above. 


The input for a sample run of BDSCAN is given in Appendix II. 


3. Restrictions Relevant to the Use of BDSCAN 


As currently written the following restrictions apply to the input data 


and problem size for BDSCAN: 


1. NAME must have < 20 characters. 
2. M must be < 60. 
3. N must be < 60. 


4. KE = 0 or J. 


To change the bounds on M and N_ you simply have to modify the dimension 
statements for A from A(61,61) to A(m' +1, n' +1) where m' and n' are 


the new bounds on M and N_ respectively. 


4. Description of Subprograms 


Detailed comments on the specific function of each subroutine may be found 
in the listing of BDSCAN in Appendix I. Variables of note are defined in 
these comments as well. In general terms: 

1. SUBROUTINE IOSUB -- Reads input cards and initiates the algorithm. 

2. SUBROUTINES MATRIX, PERMUT -- Rearrange arrays according to parameter 
instructions. 

3. SUBROUTINE FIRST1 -- Checks to see if the assumptions of the algorithms as 
specified in [1, pp. 643-644] are satisfied. 

4. SUBROUTINE PARTI -- Finds the extreme points of the region defined in 


{1, pp. 649-650] using the procedure described there. 


5. SUBROUTINE PART2 -- Finds upper bounds on the weights on these extreme points. 


{[1, pp. 650-651] 

6. SUBROUTINE PART3 -- Finds unconditional upper and lower integer bounds on 
the variables. [1l, p. 651] 

7. SUBROUTINE PART4 -- Partitions the variables into the four groups described 
in [1, pp. 651-652]. 

8. SUBROUTINE PART5 -- Obtains conditional upper and lower bounds on the 


variables in Groups 2, 3 and 4. [1, pp. 653-654] 


10. 


1A 


Ea. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


SUBROUTINE PART6 -- Initialized the procedure for the Group 1 variables. 


{l, p. 653] 

SUBROUTINE PART 7 -- Obtains a new solution for the Group 1 variables 
using the scheme described in [l, pp. 652-653]. 

SUBROUTINE PART8 -- Obtains an initial trial solution for the Group 2 
variables. [l, p. 655] 

SUBROUTINE PART9 -- Obtains a new trial solution for the Group 2 variables. 
[is ps 685) 

SUBROUTINE PART1O -- Determines if a given trial solution for the Group 2 
variables is eligible using the method described in [1l, pp. 655-658]. 
SUBROUTINE PART11 -- Initializes the bounding procedure for the variables 
in Groups 3 and 4. [1l, pp. 658-661] 

SUBROUTINE SOLVLP -- Initiates the solution of an LP via the simplex or 
dual simplex methods. These methods are described in [4]. 

SUBROUTINE SIMPLX -- Solves an LP via the simplex method. 

SUBROUTINES DUAL, INPUT, DUSEX, MP2, MP3, MP5 -- Together solve an LP via 
the dual simplex method. 

SUBROUTINE PART12 -- Obtains new lower and upper bounds for the variables 
in Groups 3 and 4. [1, pp. 660-661, 663] 

SUBROUTINE PART13 -- Obtains an initial solution for the Group 3 variables. 
{l, p. 661] 

SUBROUTINE PART14 -- Obtains a new solution for the Group 3 variables. 


{[1, pp. 661-662) 


SUBROUTINES PART15, SORT -- Together initialize the procedure for the Group 4 


variables. [1l, pp. 663-665] 


SUBROUTINES PART16 -- Obtains all the feasible completions for the Group 4 


variables. [1, pp. 663-664] 


A brief outline of how the parts of the algorithm fit together appears in 


{l, pp. 665-667). A step-by-step description of each part may be found in [2]. 


5. Description of Output 


The following output is generated by BDSCAN. During execution of SUBROUTINE 
LOSUB the following are printed. The NAME of the problem is printed out followed 
by the dimension of the constraint matrix A. Then A, the right hand side b 
and the cost coefficient array c are printed in 15F8.2 format. 

The indices of the variables that are basic in the LP solution are then 
printed followed by the LP solution itself, in 15F8.2 format. Next the suboptimal 
feasible (integer) solution that was input is printed. Finally, before execution 
of the actual algorithm is started, the number of non-slack variables that are 
basic in the LP solution is recorded [n, in l, p. 644]. 

If, during execution of SUBROUTINE FIRST1, a violation of the assumptions of 
the algorithm is found, the nature of the violation will be printed prior to 
program termination. 

On return from SUBROUTINE PART1 the first n extreme points of the region 
described in [1, p. 646] are printed in 15F8.2 format, one point at a time. 

Following SUBROUTINE PART2 the upper bounds on the weights on the first n 
extreme points are printed in 15F8.2 format. 

On return from SUBROUTINE PART3 the unconditional lower and upper bounds 
(integer) on the variables are printed in 15F8.2 format. Then the product of 
the ranges of these bounds is printed in F30.0 format. 

In SUBROUTINE PART7 a message is printed telling which component j of the 
vector X represents the variable through which eligible Group 1 solutions have 
been scanned (meaning x, has reached its upper bound and been reset to 0 in 


j 


the enumeration scheme described in [1, p. 653]). Each time j changes the 


message is printed. Also in this subroutine division by zero may result; if so, 


just prior to termination the indices K and J are printed indicating that 


XX(K,J) = 0 is the reason for the error. (Note: reference to the th entry 


in the mP extreme point printed after SUBROUTINE PART1 will confirm this.) 

During execution of SUBROUTINE SIMPLX an unbounded variable may be found. 
If so this message will be printed on return to SUBROUTINE SOLVLP, the calling 
program. 

If execution of the dual simplex method during SUBROUTINE PARTI11 finds that 
the problem being solved is infeasible, a message to that effect is printed on 
return to SUBROUTINE INPUT, the calling program. 

In SUBROUTINE PART16 each time an improved feasible (integer) solution 
is found it is printed out in 15F8.2 format, followed by its objective function 
value. 

At the end of execution of the algorithm the following information is 
printed out: 

1. The final feasibility test slacks (b - rae where xP is the optimal 
(integer) solution) in 15F8.2 format. 

2. The number of times each part was entered. 

3. The objective function value for the initial suboptimal feasible (integer) 
solution in F8.2 format. 

4. The objective function value for the LP solution in F8.2 format. 

5. The difference between the objective function value of the LP solution and 


n 
that of the ILP solution, normalized (multiplied by 1/() ae 


j=l 


J ai PLB 
format. 

6. The normalized difference between the objective function value of the ILP 
solution and that of the suboptimal integer solution in F13.5 format. 

7. The ILP solution. 


8. The objective function value of the ILP solution. 


If more than one problem was included in the input data the information 


above will be printed out for each problem starting with its input NAME. 


The output from a sample run of BDSCAN is given in Appendix III. 


6. Further Comments 

As currently written the FORTRAN code BDSCAN contains approximately 4300 
source statements. All significant array storage is done in a common block 
called COMMON. When compiled on Stanford's IBM 370/168 under IBM's FORTRAN-H- 
extended compiler, with level two optimization, approximately 512K bytes of core 
are required for all instructions and data. The code is WATFIV compatible. 


BDSCAN was last revised in November 1978. 
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LISTING OF COMPUTER PROGRAM BDSCAN 
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